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I. INTRODUCTION 



X-ray and neutron scattering techniques are probably the most effective tools to be 
employed in order to determine the structure of liquid interfaces on molecular length scales. 
These are not different in principle from conventional X-ray diffraction techniques that are 
usually applied to three dimensional crystals, liquids, solid surfaces etc. However, special 
diffractometers that enable scattering from fixed horizontal surfaces are required to carry 
out the experiments. Indeed, systematic studies of liquid surfaces had not begun until the 
introduction of the first liquid surface reflectometer. [[[[]. 

A basic property of a liquid-gas interface is the length scale over which the molecular 
density changes from the bulk value to that of the homogeneous gaseous medium. Molecular 
size and capillary waves, that depend on surface tension and gravity, are among the most 
important factors that shape the density profile across the interface and the planar correla- 
tions |I^,|9|j49[1 . In some instances the topmost layers of liquids are packed differently than 
in the bulk, giving rise to layering phenomena at the interface. Monolayers of compounds 
different than the liquid can be spread at the gas-liquid at interface, and are termed Lang- 
muir monolayers [|T7|,[5D|. The spread compound might wet the liquid surface to form a film 
of homogeneous thickness or cluster to form an inhomogeneous rough surface. The X-ray re- 
flectivity (XR) technique allows one to determine the electron density across such interfaces 
from which the molecular density and the total thickness can be extracted. The grazing 
angle diffraction (GID) technique is commonly used to determine lateral arrangements and 
correlations of the topmost layers at interfaces. GID is especially efficient in cases where 
surface crystallization of the liquid or spread monolayers occurs. Both techniques (XR and 
GID) provide structural information that is averaged over macroscopic areas in contrast to 
scanning probe microscopies (SPM's) where local arrangements are probed. For an inho- 
mogeneous interface the reflectivity is an incoherent sum of reflectivities, accompanied by 
strong diffuse scattering which in general is difficult to interpret with definitive answers and 
often requires complementary techniques to support the X-ray analysis. Therefore, prepara- 
tion of well-defined homogeneous interfaces is a key to a more definitive and straightforward 
interpretation. 



A. Competitive and Related Techniques 



Although modern scanning probe microscopies (SPM's) such as scanning tunneling mi- 
croscopy (STM) || and atomic force microscopy (AFM) || rival X-ray scattering techniques 



in probing atomic arrangements of solid surfaces, they have not yet become suitable tech- 
niques for free liquid surfaces. The large fluctuations due to the two-dimensional nature of 
the liquid interface, high molecular mobility, and the lack of electron conductivity (which is 
necessary for STM) are among some of the main obstacles that make it difficult to apply 
these techniques to gas-liquid interfaces. In dealing with volatile liquids, inadvertent depo- 
sition or wetting of the probe by the liquid can occur which may obscure the measurements. 
In addition, the relatively strong interaction of the probe with the surface might alter its 
pristine properties. For similar and other reasons, electron microscopy and electron diffrac- 
tion techniques, which are among the best choices for probing solid surfaces, are not suitable 
to most liquid surfaces, and in particular aqueous interfaces. On the other hand, visible light 
microscopy techniques such as Brewster angle microscopy (BAM) HpjUCj or fluorescence 



microscopy [KJ have been used very successfully in providing morphological pictures of the 
interface on the micrometer length scale. This information in general, is complementary to 
that extracted from X-ray scattering. These techniques are very useful for characterizing 
inhomogeneous surfaces with two or more distinct domains for which XR and GID results 
are usually difficult to interpret. However, it is impossible to determine the position of 
the domains with respect to the liquid interface, their thicknesses or their chemical nature. 
Ellipsometry is another technique that exploits visible light, to allow determination of film 
thickness on a molecular length scale and assumes that one knows the refractive index of the 
substrate and of the film |3 , 12 ; either might be different from their bulk values and difficult 
to determine. 

In the following sections theoretical background to the X-ray techniques is presented 
together with experimental procedures and data analysis concerning liquid surfaces. It is 
intended to provide a basic formulation which can be developed for further specific applica- 
tions. Several examples of these techniques applied to a variety of problems are presented 
briefly to demonstrate the strengths and limitations of the techniques. It should be borne 
in mind that the derivations and procedures described below are mostly general and can 
be applied to solid surfaces, and vice-versa, many results applicable to solid surfaces can 
be used for liquid surfaces. X-ray reflectivity from surfaces and GID have been treated in 
recent reviews H,H3|,|61I1. 



II. THEORETICAL BACKGROUND 

We assume that a plane harmonic wave of frequency uj and wave- vector k (with electric 
field, E = EQe tuJt ~ iko ' r ) is scattered from a distribution of free electrons, with a number 
density N e (r). Due to the interaction with the electric field of the X-ray wave, each free 
electron experiences a displacement proportional to the electric field, X = — -^^E. This 
displacement gives rise to a polarization P(r) distribution vector 

P(r) = N e (r)eX (1) 

in the medium. For the sake of convenience we define the scattering length density (SLD) 
p{r) in terms of the classical radius of the electron tq = 4w€ e omc 2 = 2.82 x 1CT 13 cm as follows 

p(r)=N e (r)r (2) 
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The polarization then can be written as 

The scattering length density (or the electron density) is what we wish to extract from 
reflectivity and GID experiments and relate it to atomic or molecular positions at liquid 
interfaces. The displacement vector D can now be constructed as follows: 

D = e E + P(r) = e(r)E (4) 

where e(r) is the permittivity of the medium usually associated with the refractive index 



n(r) = y e(r) . To account for absorption by the medium we introduce a phenomenological 
factor (3 that we calculate from the linear absorption coefficient fi, (given in tables plf ) as 
follows (3 = ///(2&o). Then the most general permittivity for X-rays becomes 

Air 

<r) = e Q [l-- I p{r)]+2i(3 (5) 

Typical values of the SLD (p) and the absorption term (/?) for water and liquid mercury are 
listed in Table [II In the absence of true charges in the scattering medium (i.e., a neutral 
medium) and under the assumption that the medium is nonmagnetic (magnetic permeability 
p = 1) the wave equations that need to be solved to predict the scattering from a known 
SLD can be derived from the following Maxwell equations fl3T 



V-D = V-H = 
\7xE = -dH/dt VxH = dD/dt. (6) 

Under the assumption of harmonic plane waves, E = Eoe iujt ~' lk ' r , the following general 
equations are obtained from Eq. |5|: 

V 2 E + [kl - Anp{r)]E = -V (V In e(r) • E) , (7) 
V 2 H + [kl - 4irp(r)]H = -V In e(r) x(Vxif). (8) 

In some particular cases the right hand side of Equations [7| and [S] is zero. We notice then, 
that the term 47rp(r) in the equation plays a role similar to that of a potential V(r) in 
wave mechanics. In fact, for most practical cases the right hand side of Equations [7] and |8| 
can be approximated to zero, thus the equation for each component of the fields resembles 
a stationary wave equation. In those cases general mathematical tools, such as the Born 
approximation (BA) and the distorted wave Born approximation (DWBA) can be used p? . 



A. Reflectivity 

In reflectivity experiments a monochromatic X-ray beam of wavelength A [wavevector 
k = 2tt/X and ki = (0, k y , k z )] is incident at an angle on a liquid surface and is detected 
at an outgoing angle a r such that shown in Figure [I], with a final wave vector kf. 

The momentum transfer is defined in terms of the incident and reflected beam as follows, 

Q = ki-k f , (9) 

where in the reflectivity case Q is strictly along the surface normal, with Q z = 2/c sin a = 
2k z . 
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1. Single, Ideally Sharp Interface - Fresnel Reflectivity 



Solving the scattering problem exactly for the ideally sharp interface although simple, is 
very useful for the derivation of more complicated electron density profiles across interfaces. 
The wavefunctions employed are also essential for inclusion of dynamical effects when dealing 
with non-specular scattering i.e., GID and diffuse scattering. 

la. s-polarized X-ray beam 

For a stratified medium with an electron density that varies along one direction, z, 
(p(r) = p(z) assuming no absorption i.e., (3 = 0), an s-type polarized X-ray beam with the 
electric field parallel to the surface (along the x-axis, see Fig. [I]) obeys the stationary wave 
equation as derived from Eq. |7| and is simplified as follows, 

V 2 E X + [kl - V(z)]E x = 0, (10) 

with an effective potential V(z) = 47cp(z). The general solution to Eq. is then given by, 

E x = E(z)e ik * y (11) 

where the momentum transfer along y is conserved when the wave travels through the 
medium leading to the well-known Snell's rule for refraction. Inserting Eq. |TT| in Eq. |It] 
leads to a one dimensional wave equation through a potential V(z), 

+ \k\ - V(z)}E = 0. (12) 

The simplest case of Eq. [12| is that of an ideally flat interface, separating the vapor phase 
and the bulk scattering length density p s , at z = 0. The solution of Eq. is then given by 

F( \ = I eikz '° Z + r ( k z,s)e~ ihzfiZ z > in gas , . 

[Z) \ t(k ZtS )e ik ^ z z<0m liquid 1 j 

where 



Ks = V k lo ~ 47r Ps = V k lo - k l (14) 

where k c = 1^J-np s . By applying continuity conditions to the wavefunctions and to their 
derivatives at z — 0, the Fresnel equations for reflectance, r(fc ZjS ), and transmission, t(k z>s ), 
are obtained, 

r(M = *(M = (15) 

The measured reflectivity from an ideally flat interface, Ftp, is usually displayed as a function 
of the momentum transfer Q z = k z $ + k z )S pa 2k z , and is given by 

R F (Q z ) = \r(k z>s )\ 2 . (16) 

Below a critical momentum transfer, Q c = 2k c = 4w^7fpJ, k z ^ s is an imaginary number and 
Rf{Qz) — 1; total external reflection occurs. Notice that whereas the critical momentum 
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transfer does not depend on the X-ray wavelength, the critical angle for total reflection does, 
and it is given by a c ~ xJ p s /ir. Typical values for critical angles for X-rays of wavelength 

^CuK a = 1.5404A, are listed in Table |. For Q z ^> Q c , Rf(Q z ) can be approximated to a 
form that is known as the Born approximation, 

This form of the reflectivity at large Q z s is also valid for internal scattering, i.e., reflec- 
tivity from liquid into the vapor phase. However, total reflection does not occur for the 
internal reflectivity case. Calculated external and internal reflectivity curves from an ideally 
flat surface, Rf, displayed versus momentum transfer (in units of the critical momentum 
transfer) are shown in Fig. 0(a). Both reflectivities converge at large momentum transfer 



where they can be both approximated by Eq. 17. The dashed line in the same figure shows 



the approximation (Q C /2Q Z ) A , which fails in describing the reflectivity close to the critical 
momentum transfer. 

The photon transmission at a given k zs is given by 



|2 Re(fc 2 



nh,) = \tiK,)\'^y (18) 

where the ratio on the right hand side accounts for the flux through the sample. In the case 
of external reflection, and for values of k z ^ that are smaller than k c , the real part of k z>s 
is zero, and there is no transmission, whereas above the critical angle k z ^ s is real, and the 
transmission is given by, 

T{k z , s ) = ^'° fc ;% for k z , > k c (19) 
{k zfi + k z |- 



1 Z,SJ 



and the conservation of photons is fulfilled in the scattering process, 

T{k z>s ) + R(k Z)S ) = 1. (20) 

In Fig. H](b) the transmission amplitude |£(& Z)S )| for external (solid line) and for internal 
(dashed line) reflections are shown. This amplitude modulates non-specular scattering pro- 
cesses at the interface as will be discussed later in this unit. 

The effect of absorption on the reflectivity can be incorporated by introducing f3 into 



the generalized potential in Eq. |TTj, so that k z>s = J k z — k% + 2i(3 is used in the Fresnel 
equations Eq. |15[ Calculated reflectivities from water and liquid mercury demonstrating 
that the effect of absorption is practically insignificant for the former yet has the strongest 
influence near the critical angle for the latter, are shown in Fig. 0. 



2. Multiple Stepwise and Continuous Interfaces 

On average, the electron density across a liquid interface is a continuously varying func- 
tion, and is a constant far away on both sides of the interface, as is shown in Fig. 01. 
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The reflectivity for a general function p{z) can be then calculated by one of several 
methods classified into two major categories: dynamical and kinematical solutions. The 
dynamical solutions are in general more exact and include all the features of the scattering, in 
particular the low angle regime, close to the critical angle where multiple scattering processes 
occur. For a finite number of discrete interfaces exact solutions can be obtained by use of 
standard recursive |36| or matrix |7|] methods. These methods can be extended to compute 



with very high accuracy the scattering from any continuous potential by slicing it into a 
finite but with sufficient number of interfaces. On the other hand, the kinematical approach 
neglects multiple scattering effects and fails in describing the scattering at small angles. 

2a. The Matrix Method 
In this approach the scattering length density with variation over a characteristic length 
dt is sliced into a histogram with N interfaces. The matrix method, is practically equivalent 
to the Parratt formalism 0. For each interface, the procedure described previously for the 
one interface is applied. Consider an arbitrary interface, n, separating two regions, of a sliced 
SLD (as in Fig. f|), with p n _i, and p n at position z = z n with the following wavefunctions 



Rn 



2&n — 1 Z 



Pn-l 

i^n— 1 Z 



Pn 

D -ik n z 

I: A - -f f Z 



rri 



(21) 



where L, 



k 2 zQ — Anp n . The effect of absorption can be taken into account as described 



earlier. For simplicity, the subscript z is omitted from the component of the wavevector so 
that k zn = k n . The solution at each interface in terms of a transfer matrix, M n , is given by 




ip ^{kn — 1 ~\~k n ^jz n ^i{fcn— 1 k n ^jz n 




where r n = 1 ^" is the Fresnel reflection function through the z n interface separating the 
p n -i and p n SLD's. The solution to the scattering problem is given by noting that beyond 
the last interface, (i.e., in the bulk), there is a transmitted wave only for which an arbitrary 
amplitude of the form (J) can be assumed (i.e., the reflectivity is normalized to the incident 
beam anyway). The effect of all interfaces is calculated as follows 

(SaH M ') (M2) "■(""«) (J) 

with rjy + i = in the Mn+i matrix given in terms of the substrate k s . The reflectivity 

is then given by the ratio 



R(Qz 



Rq,i 



(22) 



Applying this procedure to the one box model of thickness d with two interfaces yields 



R{Q Z = 2k s 



n + r 2 e 
1 + rir 2 e 



i2k s d 



i2k s d 



(23) 
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Figure |5| shows the calculated reflectivities from a flat liquid interface with two kinds of 
films (one box) of the same thickness d but with different scattering length densities, p\ and 
Pi- The reflectivities are almost indistinguishable when the normalized SLD's {pi/p s ) of the 
films are complementary to one (p\/p s + P2/ Ps — 1), except for a very minute difference near 
the first minimum. In the kinematical method described below, the two potentials shown in 
Fig. [5] yield identical reflectivities. 

The matrix method can be used to calculate the exact solution from a finite number 
of interfaces, and it is most powerful when used with computers by slicing any continuous 
scattering length density into a histogram. The criteria for determining the optimum number 
of slices to use is based on the convergence of the calculated reflectivity at a point where 
slicing the SLD into more boxes does not change the calculated reflectivity significantly. 

2b. The Kinematical Approach 

The kinematical approach for calculating the reflectivity is only applicable under certain 
conditions where multiple scattering is not important. It usually fails in calculating the 
reflectivity at very small angles (or small momentum transfers) near the critical angle. The 
kinematical approach, also known as the Born- Approximation, gives physical insight in the 
formulation of R(Q Z ) by relating the Fresnel normalized reflectivity, R/Rp, to the Fourier 
transform of spatial changes in p(z) across the interface as discussed below. 

As in the dynamical approach, p(z) is sliced so that k(z) = 
across an arbitrary point z is given by 



^2,0 



k 2 and the reflectance 



r(z) 



k(z + Az) - k(z) 



4n dp 



k(z + Az) + k(z) 4k(z) 2 dz 



dz « (Q c /2Q z ) 2 -^dz 
p s dz 



(24) 



In the last step of the derivation, r(z) was multiplied and divided by p s , the SLD of the 
subphase, and the identity Q 2 = I6np s was used. Assuming no multiple scattering, the 
reflectivity is calculated by integrating over all reflectances at each point, z, with a phase 
factor e lQzZ as follows, 



R(Qz) = Rf(Qz 



Rf(Qz)\®(Qz 



(25) 



where <&(Q Z ) can be regarded as the generalized structure factor of the interface, analogous 
to the structure factor of a unit cell in 3D crystals. This formula can be also derived by 
using the Born Approximation, as is shown in the following section. 

As an example of the use of Eq. |25| we assume that the SLD at a liquid interface can be 
approximated by a sum of error functions as follows 



P(z) = Po + 2^ 

3=1 Z 



1 + erf 



( z ~ z i 
\V2ai 



(26) 



where po is the SLD of the vapor phase and pn = p s - Using Eq. |25| the reflectivity is given 
by 



R(Qz) = Rf(Qz) 



E 



- Ps . 



exp 



exp 



(27) 
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Assuming one interface at Z\ = with surface roughness o\ = a, the Fresnel reflectivity, 
Rp(Q z ), is simply modified by a Debye- Waller-like factor 



R(Q z ) = R F (Q z )exp-^ 2 . (28) 

The effect of surface roughness on the reflectivities from water and from liquid mercury 
surfaces assuming gaussian smearing of the interfaces, is shown in Figure ||. Braslau et al. 
|| have demonstrated that the Gaussian smearing of the interface due to capillary waves in 
simple liquids is sufficient in modeling the data, and that more complicated models cannot 
be supported by the X-ray data. 

Applying Eq. |27] to the one box model discussed above (See Fig. |]), and assuming 
conformal roughness, o~j = a, the calculated reflectivity in terms of SLD normalized to p s is 

R{Qz) = Rf(Qz) [(1 - Pi) 2 + P\ + 2pi(l - Pi) cos(Q z d)] e^^\ (29) 

In this approximation the roles of the normalized SLD of the one box, p 1; and that for the 
complementary model P2 = 1 — p\ are equivalent and demonstrates that the reflectivities 
for both models are mathematically identical. This is the simplest of many examples where 
two or more distinct SLD models yield identical reflectivities in the Born Approximation. 
When using the kinematical approximation to invert the reflectivity to SLD there is always a 
problem of facing a non- unique result. Ways to distinguish between such models is discussed 
in the Data Analysis section. 

In some instances, the scattering length density can be generated by several step functions 
that are smeared with one gaussian {conformal roughness Oj = a), representing different 
moieties of the molecules on the surface. The reflectivity can be calculated by using a 
combination of the dynamical and the kinematical approaches ||. First, the exact reflectivity 
from the step-like functions (a = 0) is calculated using the matrix method, Rd yn (Qz), and 
the effect of surface roughness is incorporated by multiplying the calculated reflectivity with 
a Debye- Waller-like factor as follows f|| 

R(Q z ) = R dyn (Q z )exp^^ 2 . (30) 



B. Non-specular scattering 

The geometry for non-specular reflection is shown in Figure [L](b) . The scattering from 
a 2D system is very weak and enhancements due to multiple scattering processes at the 
interface is taken advantage of. As is shown in Figure |l|(b) the momentum transfer Q has 
a finite component parallel to the liquid surface (Q± = k\ — k ± ) enabling determination 
of lateral correlations in the 2D plane. Exact calculation of scattering from surfaces is 
practically impossible except for special cases, and the Born Approximation (BA) |47|] is 
usually applied. When the incident beam or the scattered beam are at grazing angles 
(i.e., near the critical angle), multiple scattering effects modify the scattering and these 
can be accounted for by a higher order approximation known as the Distorted Wave Born 
Approximation (DWBA). The features due to multiple scattering at grazing angles provide 
evidence that the scattering processes indeed occur at the interface. 
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1. The Born Approximation 



In the BA for a general potential V(r) the scattering length amplitude is calculated as 



follows m 



F(Q) = -!- J V(r)e iQ - r d 3 r. 



(31) 



where in the present case, V(r) = Airp{r). From the scattering length amplitude the differ- 
ential cross section is calculated as follows 



[p{0)p{r)] e iQ - r d 3 r. 



(32) 



where [p(0)p(r)] = f[p(r' — r)p(r')]d 3 r' is the density-density correlation function. The 
measured reflectivity is a convolution of the differential cross-section with the instrumental 
resolution, as discussed below and in the literature [HTL|S|,ff9 . 



The scattering length density, p, for a liquid-gas interface can be described as a function 
of the actual height of the surface, z(x,y), as follows, 



p{fi,z) 



p s for z < z(/j,) 
for z > z(fjt) 



(33) 



where [i = (x, y, 0) is a 2D in-plane vector. The height of the interface z is also time depen- 
dent and temperature dependent due to capillary waves, and therefore thermal averages of 



z are used [|10| , |T4|| . Inserting the SLD Eq. |33], in Eq. ^ and performing the integration over 
the z coordinate yields 



F(Q±,Q Z ) = / e^^^V/x 



(34) 



where Q± = (Q x , Qy, 0) is an in-plane scattering vector. This formula properly predicts the 
reflectivity from an ideally flat surface, z(x, y) = within the kinematical approximation 



F(Qx,Q z ) 



4vr 2 p. 



iQz 



(35) 



with a 2D 6 function that guarantees specular reflectivity only. The differential cross-section 
is then given by 



da 
dfi 



7T 



'01 
AQz 



(36) 



where Q 2 C = 16irp s . This is the general form for the Fresnel reflectivity in terms of the 
differential cross-section da/dQ, which is defined in terms of the flux of the incident beam 
on the surface. In reflectivity measurements however, the scattered intensity is normalized 
to the intensity of the incident beam and therefore the flux on the sample is angle dependent 
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and is proportional to sinctj. In addition the scattered intensity is integrated over the polar 
angles af and 29 with k$ sin a i dajd(26 l ) = dQ x dQ y , correcting for the flux and integrating 



R F (Qz) « (TtatiQz) = JJ[-£ 



da\ dQxdQy _ ( Q c 
dVt ) Air 2 kQ sin ai sin a f \2Q Z 



(37) 



as approximated from the exact solution, given in Eq. |TT| . 

Taking advantage of the geometrical considerations above, the differential cross-section 
to the reflectivity measurement can be readily derived in the more general case of scattering 
length density that varies along z only, (i.e., p(z)). In this case, Eqs. |31] can be written as 



da 

dn 



47rV(Q, 



p{z)e iQ * z dz 



(38) 



If we normalize p(z) to the scattering length density of the substrate, p s , and use a standard 
identity between the Fourier transform of a function and its derivative, we obtain, 



da 

dn 



71 



'01 

4Q, 



1 dp(z) 
Ps dz 



*dz 



(39) 



which with the geometrical corrections yields Eq. |25|. 

Thermal averages of the scattering length density under the influence of capillary waves 
and the assumption that the SLD of the gas phase is zero can be approximated as follows 



[p(0)p(r)] 



P ~ l + erf[^— 
2 V [ V2a(fi) 



(40) 



where cr(/z) is the height-height correlation function. 

Inserting Eq. HO into Eq. |32| and integrating over z, results in the differential cross-section 



da 

dh~ 



2tt 



jQx -fj,- 



(41) 



and assuming isotropic correlation function in the plane yields ||9] 

pJ (Q±^>- QWM dp 



da 

dh~ 



2tt 



(42) 



where Jo is a Bessel function of the first kind. This expression was used by Sinha et al. 
to calculate the diffuse scattering from rough liquid surfaces with a height-height density 
correlation function that diverges logarithmically due to capillary waves [fE3.[It] . 



2. Distorted Wave Born- Approximation (DWBA) 

Due to the weak interaction of the electromagnetic field (X-rays) with matter (electrons), 
the BA is a sufficient approach to the scattering from most surfaces. However, as we have 
already encountered with the reflectivity, the BA fails (or is invalid) when either the incident 
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beam or the scattered beam is near the critical angle, where multiple scattering processes 
take place. The effect of the bulk on the scattering from the surface can be accounted for 
by defining the scattering length density as a superposition of two parts as follows 

p(r)= Pl (z)+p 2 (v,z). (43) 

Here p\ (z) is a step function that defines an ideally sharp interface separating the liquid and 
gas phases at z — 0, whereas the second term, p 2 (fi, z), is a quasi two-dimensional function in 
the sense that it has a characteristic average thickness d c such that lim z _ > j-d e /2 P2(a*> z ) = 0- 
It can be thought of as film-like and is a detailed function with features that relate to 
molecular or atomic distributions at the interface. Although the definition of p 2 may depend 
on the location of the interface, (z = 0) in pi, the resulting calculated scattering must be 
invariant for equivalent descriptions of p{r). In some cases p 2 can be defined as either a 
totally external or totally internal function with respect to the liquid bulk, (i.e., pi). In 
other cases, especially when dealing with liquid surfaces, it is more convenient to locate the 
interface at some intermediate point coinciding with the center of mass of p 2 with respect to 
z. The effect of the substrate term P i (z) on the scattering from p 2 can be treated within the 
distorted wave Born approximation (DWBA) by using the exact solution from the ideally 
flat interface (Sec. II) to generate the Green function for a higher order Born approximation 
P?| , |4^|J57| , |4^[] . The Green function in the presence of an ideally flat interface, pi, replaces 
the free particle Green function that is commonly used in the Born approximation. The 
scattering amplitude in this case is given by [47.42f 



Fdwba(Q) = F F (Q Z ) + F 2 (Q) = -inQ z r F (Q z ) + J xUr)p2{r) X h{r)dr (44) 

where the exact Fresnel amplitude Fp{Q z ) is written in the form of a scattering amplitude 
so that Eq. ^4| reproduces the Fresnel reflectivity in the absence of p 2 . 
The exact solution of the step function pi, Xk( r ), is given by 



Xfe(r) " e < for^>0 



and the Xk( r ) * s the time reversed and complex conjugate solution of an incident beam with 
-*/ \ -ik f -u I t f (k f z c)e~ lk ^° z for z > 

Xt(r) — e -L ^ J v *fl' (46) 
| e-<° z + rf{k f Zt0 )e<°* for z < 0. 

In the absence of pi, Eq. P| transforms into the standard Born Approximation for p 2 and this 
reproduces Eq. [31]. The notation for transmission and reflection functions indicate scattering 
of the wave from the air onto the subphase and vice-versa according to k l z s = J {k l z ) 2 — k^, 

and k[ = \j (kz,s) 2 + k% respectively. In the latter case, total reflectivity does not occur 
except for the trivial case k{ s = 0, and no enhancement due to the evanescent wave is 
expected. In this approximation, the final momentum transfer in the Q z direction is a 
superposition of momentum transfers from p 2 (the film) and from the liquid interface, p\. 
For instance, there could be a wave scattered with Q z = with respect to p 2 but reflected 
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from the surface with a finite Q z . Assuming that the scattering from the film is as strong 
for Q z as for — Q z (as is the case for an ideal 2D system with equal scattering along the rod, 
i.e., p 2 is symmetrical under the inversion of z), we can write p~2(Q±,Qz) ~ p 2 (Qj_,— Q z ). 
It can be shown that the cross-section can be approximated as follows |57|JT5|j49[1 



da 



m tS )p 2 (Q ± ,Q' z ) t f(ki s ) 



2 

for z > 

(47) 

for z < 



where Q± = k l ± — kj_ and Q' z = k l zQ — k{ s and Q" = k\ s — k( . Notice that the transmission 
functions modulate the scattering from the film (P2), and in particular they give rise to 
enhancements as k\ s and k{ s are scanned around the critical angle as depicted in Figure 
0(b). Also, it is only by virtue of the z symmetry of the scatterer that such enhancements 
occur for an exterior film. From this analysis we notice that there will be no enhancement 
due to the transmission function of the final wave for an interior film. 

To examine the results from the DWBA method, we consider scattering from a single 
scatterer near the surface. The discussion is restricted to that where the detection of the 
scattered beam is performed in the vapor phase only. The scatterer can be placed either 
in the vacuum (z > 0) or in the liquid (see Fig. ||). When the particle is placed in the 
vacuum there are two major relevant incident waves: a direct one from the source, labeled 
1* in Figure |6|(a); and a second one reflected from the surface before scattering from p 2 , 
labeled 2*. Assuming inversion symmetry along z, both waves scatter into a finite Q± with 
similar strengths at Q z and — Q z , giving rise to an enhancement near the critical angle 
if the incident beam is near the critical angle. Another multiple scattering process that 
gives rise to enhancement at the critical angle is one in which beam 1* does not undergo a 
change in the momentum transfer along (z Q{ dm ~ 0) before scattering from the liquid 
interface. The effect of these processes gives rise to enhancements if either the incident 
beam or the reflected beam are scanned along the z direction. Slight modifications of the 



momentum transfer along the z direction (such as Q' z = k\ + J \kl) 2 — k 2 ) are neglected 
in the discussion above). The effective amplitude from the scatterer outside the medium is 
given by the following terms 

[e «?.* + e «*Mfc£j] « t(K, s ) (48) 

where the approximation is valid since at small momentum transfers the phase factor can 
be neglected, and 1 + r(k Z)S ) = t(k z>s ). At large angles the reflectivity is negligible and 
the transmission function approaches t(k ZyS ) ~ 1. Similar arguments hold for the outgoing 
wave. Neglecting the small changes in the momentum transfer due to dynamical effects, the 
transmission function modulates the scattering as is shown in solid line in Fig. @(b). 

The scattered wave from a particle that is embedded in the medium is different due 
to the asymmetry between external and internal reflection from the liquid subphase. The 
wave scattered from the particle re-scatters from the liquid-gas interface. Upon traversing 
the liquid interface the index of refraction increases from n = 1 — 5 (5 = (27t//cq)p to 1 
and no total internal reflection occurs as discussed earlier; thus there is no evanescent wave 
in the medium. The transmission function for this wave is given by t(—k l z J, like that of a 
wave emanating from the liquid interface into the vapor phase. In this case, the transmission 
function is a real function for all k l z s and does not have the enhancements around the critical 
angle (as shown in Figure §(b)), with zero intensity at the horizon. 
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3. Grazing Incidence Diffraction ( GID ), and Rod Scans 



In some instances ordering of molecules at liquid interfaces occurs. Langmuir monolayers 
spread at the gas-water interface usually order homogeneously at high enough lateral pres- 
sures. |^,|n|^4]]. Surface crystallization of n-alkane molecules on molten alkane has been 
observed recently |p5| , pi| . In these cases, p 2 is a periodic function in x and y, and can be 
expanded as a Fourier series in terms of the 2D reciprocal lattice vectors t± as follows 

p 2 (LL,z)=Y,F(r ± ,z)e i ^ (49) 



Inserting Eq. ^ in Eq. 47 and integrating yields the cross-section for quasi-2D Bragg reflec- 
tion at = t ± 

^ ~ P(Q) \t(kl s )\ 2 (\F(t ± , Q z )\ 2 ) DW(Q ± , Q z ) \t f (k{ s )\ 2 8{Q - tjl) (50) 

where P{Q) is a polarization correction, and the 2D unit cell structure factor is given as a 
sum over the atomic form factors fj(Q) with appropriate phase 

F(T ± ,Q z ) = Y,fAQ> lT - rj+QzZj - (51) 

3 

The structure factor squared is averaged for multiplicity due to domains and weighted for 
orientation relative to the surface normal. The ordering of monolayers at the air-water 
interface is usually in the form of 2D powder consisting of crystals with random orientation 
in the plane. From Eq. |5D] we notice that the conservation of momentum expressed with 
the 8 function allows for observation of the Bragg reflection at any Q z . A rod scan can be 
performed by varying either the incident or reflected beam, or both. The variation of each 
will produce some modulation due to both the transmission functions and to the average 
molecular structure factor along the z-axis. The Debye- Waller factor , DW(Q±, Q z ), which 
is due to the vibration of molecules about their own equilibrium position with time dependent 
molecular displacement u(t) is given by 

DW(Q X , Qz) ~ e -^Ql<^l>+Q> 2 ) (52) 

The term due to capillary waves on the liquid surface is much more dominant than the 
contribution from the in-plane intrinsic fluctuations. The Debye- Waller factor in this case 
is an average over a crystalline size and might not reflect surface roughness extracted from 
reflectivity measurements, where it is averaged over the whole sample. 



III. EXPERIMENTAL CONSIDERATIONS AND DATA ANALYSIS 



The minute sizes of interfacial samples on the sub-microgram level combined with the 
weak interaction of X-rays with matter result in very weak GID and reflectivity (at large Q z ) 
signals that require highly intense incident beams, which are available at X-ray synchrotron 
sources. A well prepared incident beam for reflectivity experiments at a synchrotron (for 
example, the X22B beam-line at the National Synchrotron Light Source at Brookhaven 
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National Laboratory) has an intensity of 10 9 — 10 10 photons/second, whereas for a similar 
resolution, an 18 kW rotating anode generator produces 10 4 — 10 5 photons/second. Although 
reflectivity measurements can be carried out with standard X-ray generators, the measure- 
ments are limited to almost half the angular range accessible at synchrotron sources and they 
take hours to complete compared to minutes at the synchrotron. GID experiments are prac- 
tically impossible with X-ray generators since the expected signals (2D-Bragg reflections, 
for example) normalized to the incident beam are on the order of 10~ 8 — 10 -10 . 



A. Reflectivity 

X-ray reflectivity and GID measurements of liquid surfaces are carried out on special 
reflectometers that enable manipulation of the incident as well as the outgoing beam. A 
prototype liquid surface reflectometer was introduced for the first time by Als-Nielsen and 
Pershan |TJ. In order to bring the beam to an angle of incidence «j with respect to the liquid 
surface, the monochromator is tilted by an angle x either about the axis of the incident beam 
(indicated by xi i n Fig. 0) or about the axis normal to the reciprocal lattice wave vector of 
the monochromator, r {xi)- Figure |7| shows the geometry that is used to deflect the beam 
from the horizon onto the liquid surface at an angle «j by tilting the monochromator. At the 
Bragg condition, the surface of the monochromator crystal, is at an angle tp with respect to 
the incoming beam. Tilting over the incident beam axis is like tracing the Bragg reflection 
on the Debye-Scherer cone so that the ip axis remains fixed, with a constant wavelength at 
different tilting angles. The rotated reciprocal lattice vector and the final wave-vector in 
this frame are given by 

T o = To( — si n i>, cosip cosxi, cosip sin xi) 

kf = /c (cosa;i cos0, cosctj sin0, sinaj) (53) 
where <fi is the horizontal scattering angle. The Bragg conditions for scattering are given by, 

ki + T = kf, \k f \ = k (54) 
Using Eqs. [53] and 153 the following relations for the monochromator axes are obtained, 



• / r ° 
sin w = —-- 

2k 

sinxi = — cos?/' sin «j (55) 
to 

cos0 = (l - ^\ I cos a* 



And we notice that the monochromator angle ip is independent of c^. However, the scattering 
angle <fr has to be modified as is varied. This means that the whole reflectometer arm has 
to be rotated. Similarly, for the configuration where the monochromator is tilted over the 
axis normal to r we get 
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sm ^ = 7Tj 

2k cos X2 

^0 

sinY 2 = — sinaj (56) 

cos0 = ^1 - I cos a». 

From these relations the conditions for a constant wavelength operation for any angle of 
incidence, a i} can be calculated and applied to the reflectometer. Here, unlike the previous 
mode, deflection of the beam to different angles of incidence requires both, the adjustment 
of if) as well as in order to maintain a constant wavelength. If if) is not corrected in this 
mode of operation, the wavelength varies as x is varied. This mode is sometimes desir- 
able, especially when the incident beam hitting the monochromator consists of a continuous 
distribution of wavelengths around the wavelength at horizontal scattering, \2 — 0. Such 
continuous wavelength distribution exists when operating with X-ray tubes, or when the 
tilting monochromator is facing the white beam of a synchrotron. Although, the variation in 
the wavelength is negligible as \2 is varied, without the correction of if>, the exact wavelength 
and the momentum transfer can be computed using the relations in Eq. [5^]. In both modes 
of monochromator tilting, the surface height as well as the height of the slits are adjusted 
with vertical translations. 

The exact angle of the monochromatic incident beam on the surface is determined by 
at least two horizontal slits located between the sample and the source. One of these slits 
is usually located as close as possible to the sample, and the other as close as possible to 
the source. These two slits determine the resolution of the incident beam. By deflecting the 
beam from the horizon the shape of the beam changes and that may change incident beam 
intensity going through the slits, and the use of a monitor right after the slit in front of the 
sample is essential for the absolute determination of the reflectivity. The size of the two slits 
defining the incident beam, is chosen in such a way that the foot-print of the beam is much 
smaller than the width of the reflecting surface so that total reflectivity occurs. Figure |8| 
shows the reflected beam and the direct beam from a flat surface of water demonstrating 
total reflectivity at Q z = 0.85<5 C - In this experiment the detector slit is wide open at about 
ten times the opening of the sample slit. As is demonstrated, the effect of absorption is 
negligible for water, and roughness is significantly reduced by damping surface waves. The 
damping can be achieved by reducing the height of the water film to about ~ 0.3mm and 
placing a passive as well as an active anti- vibration unit underneath the liquid sample holder, 



suppressing mechanical vibrations p4 



B. Non-specular Scattering - GID, Diffuse Scattering and Rod Scans 

X-ray GID measurements are performed at angles of incidence below the critical angle 
~ 0.9a c . Operating with the incident beam below the critical angle enhances the signal from 
the surface with respect to that of the bulk, by creating an evanescent wave in the medium 
that is exponentially decaying as, 

E{z) = t(k z , s )e~ z / A (57) 
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where 




For water at k z « 0.9k c , A ~ 1000 A. 

As illustrated in Figure 0(b), the components of the momentum transfer for GID are 
given by, 

Qz = k (smai + sin af) 

Qx = /co(cos ati — cos a/) cos 26* 

Q y = k cos a f sin 29 (59) 

In most cases, the 2D order on liquid surfaces is powder like, and the lateral scans are 
displayed in terms of Q± which is given by, 

Q± = k \J cos 2 OLi + cos 2 OLj — 2 cos a j cos a f cos 29 (60) 

To determine the in-plane correlations the horizontal resolution of the diffractometer can 
be adjusted with a soller collimator that consists of vertical absorbing foils stacked together 
between the surface and the detector. The area that is probed at each scattering angle 29 is 
proportional to Sq/ sin 29, where So is the area probed at 29 = tt/2. The probed area must 
be taken into account in the analysis of a GID scan that is performed over a wide range of 
angles. 

Position sensitive detectors (PSD) are commonly used to measure the intensity along the 
2D rods. It should be kept in mind that the intensity along the PSD is not a true rod scan 
of a Bragg reflection at a nominal Q± because of the variation in Q± as a/ is varied as is 
seen in Eq. |59|. 



C. Data Analysis 

The task of finding the SLD from a reflectivity curve is similar to that of finding an 
effective potential for Eq. |7| from the modulus of the wave-function. Direct inversion of the 
scattering amplitude to SLD is not possible except for special cases when the BA is valid [03 



If the modulus and the phase are known they can be converted by the method of Gelfand- 
Levitan-Marchenko [jPJ] to SLD (GLM method). However, in reflectivity experiments the 
intensity of the scattered beam alone is measured, and phase information is lost. 

Direct reconstruction of step-like potentials have been developed recently by retrieving 
the phase from the modulus i.e., reflectivity and then using the GLM method fflJTT]. Model- 
independent methods which are based on optimization of a model to reflectivity, without 
the requirement of any knowledge of the chemical composition of the SLD at the interface 
were also developed recently [p?| , |6T| . Such models incorporate a certain degree of objectivity. 



These methods are based on the kinematical and the dynamical approaches for calculating 



the reflectivity. One method |37j uses indirect Fourier transformation to calculate the cor- 



relation function of dp/dz which is subsequently used in a square-root deconvolution model 
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to construct the SLD model. Zhou and Chen, on the other hand, developed a groove track- 
ing method that is based on an optimization algorithm to reconstruct the SLD using the 
dynamical approach to calculate the reflectivity at each step |61 . 



The most common procedure to extract structural information from reflectivity is by 
using standard non-linear least squares refinement of an initial SLD model. The initial 
model is defined in terms of a P-dimensional set of independent parameters, p, using all the 
available information in guessing of p(z, p). The parameters are then refined by calculating 
the reflectivity (R[Q l z ,p\) with the tools described earlier and by minimizing the x 2 {p) 
quantity, 



R exp (Ql)-R(Ql,p) 
<Ql) 



(61) 



where e{Q\) is the uncertainty of the measured reflectivity, R exp {Q l z ) and N is the number 
of measured points. The criteria for a good fit can be found in [Q. Uncertainties of a certain 
parameter can be obtained by fixing it at various values and for each value refining the rest 
of the parameters until \ 2 is increased by a factor of at least 1/(N — P). 

The direct methods and model-independent procedures of reconstruction SLD do not 
guarantee uniqueness of the potential i.e., there can be multiple SLD profiles that essentially 
yield the same reflectivity curve, as discussed with regard to Figure |5], for example. The 
uniqueness can be achieved by introducing physical constraints that are incorporated into 
the parameters of the model. Volume, in-plane density of electrons etc., are among such 
constraints that can be used, (applying such constraints is discussed briefly in the Examples 



section |51|j5^ , p!8[| ) . These constraints reduce the uncertainties and make the relationship of 
the SLD to the actual molecular arrangement apparent. In the dynamical approach no two 
potentials yield exactly the same reflectivity although the differences between two models 
might be too small to be detected in an experiment. 

An experimental method to solving such a problem was suggested by Sanyal et al., using 
anomalous X-ray reflectivity methods. Two reflectivity curves from the same sample are 
measured with two different X-ray energies, one below and one above an absorption edge 



of the substrate atoms, thereby varying the scattering length density of the substrate |}46 
Subsequently the two reflectivity curves can be used to perform a direct Fourier reconstruc- 
tion j5| or by refinement methods to remove ambiguities. This method is not efficient when 
dealing with liquids that consist of light atoms because of the very low energy of the absorp- 
tion edge with respect to standard X-ray energies. Another way to overcome the problem 
of uniqueness is by performing reflectivity experiments on similar samples with X-rays and 
with neutrons. In addition, the SLD, p(z) across the interface can be changed significantly, 
in neutron scattering experiments, by chemical exchange of isotopes that change p(z), but 



maintain the same structure ||51|| . The reflectivities (X-ray as well as neutrons) can be fitted 
to one structural model that is defined in terms of geometrical parameters only, calculating 
the SLD's from scattering lengths of the constituents and the geometrical parameters pT|.|^ 



IV. EXAMPLES 

Since the pioneering work of Als-Nielsen and Pershan QXJ] , X-ray reflectivity and GID be- 
came standard tools for the characterization of liquid surfaces on the atomic length scales. 
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The techniques have been exploited in studies of the physical properties of simple liquids 
Lanemuir monolayers g|]|[|||4],|,|]] , liquid metals |l]|n],|4| 
liquid crystals 



surface crystal- 
surface properties of quantum liquids pj, protein 



lization 

recognition processes at liquid surfaces |53],|54|,|29|] and many other applications. Here, only 
several examples are briefly described in order to demonstrate the strengths and the limi- 
tations of the techniques. In presenting the examples, there is no intention of giving a full 
theoretical background of the systems. 



A. Simple Liquids 

The term simple liquid is usually used for mono-atomic systems governed by Van der 
Waals type interactions such as, liquid argon. Here, the term is extended to include all clas- 
sical dielectric (non-metallic) liquids such as water, solvents (methanol,ethanol, chloroform 
etc.) and others. One of the main issues regarding dielectric liquids is the determination of 
the average density profile across the interface, N(z). This density is the result of folding the 
intrinsic density iVj(z) of the interface due molecular size, viscosity, and compressibility of 
the fluid with density fluctuations due to capillary waves, 6Ncw{z)- The continuous nature 
of the density across the interface due to capillary waves was worked out by Buff, Lovett and 
Stillinger (BLS) |T(J assuming that Nj(z) is an ideal step like function. The probability for 
the displacement is taken be to proportional to the Boltzmann factor, e~^ u< ^ z \ where U is 
the free energy necessary to disturb the surface from equilibrium state (i.e., z(x, y) =0), and 
(3 = 1/ksT. The free energy of an incompressible and non- viscous liquid surface consists of 
two terms; a surface tension (7) term, that is proportional to the changes in area from the 
ideally flat surface and a gravitational term as follows, 

U = f(i yJl + \Vz\ 2 -l +\m s gz 2 ^& 2 ii 

|/( 7 |Vz| 2 + m s ^ 2 )dV (62) 



~ 2- 



where m s is the mass density of the liquid substrate. By using standard Gaussian approxi- 
mation methods Buff et al find, that U(z) ~ ^-7. After convolution of the probability with a 
step like function, representing the intrinsic density of the liquid surface, yields the following 
density function, 

N(z) = A s erfc (-^j (63) 

with a form similar to the one given in Eq. The average surface roughness at temperature 
T, is then given by 

where is a molecular diameter and L is the size of the surface. Notice the logarithmic 
divergence of the fluctuations as the size of the surface increases, as expected of a 2D system 



|25[. This model was further refined by assuming that the intrinsic profile has a finite width 
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14]. In particular if the width due to the intrinsic profile is also expressed by a Gaussian 



then, the effective surface roughness is given by 

°eff = ° 2 i + °lw ( 65 ) 

and the calculated reflectivity is similar to Eq. |28| for an interface that is smeared like the 
error function. 

Row = R F (Q z )e-<ff Q * (66) 

Fig. ^] shows the reflectivity from pure water measured at the synchrotron BH] where it 
is shown that using Eq. ^6] for fitting the reflectivity data is satisfactory implying that the 
error function type of density profile (BLS model) for the liquid interface is sufficient. In the 
refinement procedure only one parameter, the surface roughness a, is varied (cr = 2.54 A). 
This small roughness value depends on the attenuation of capillary waves by minimizing the 
depth of the water to about 0.3mm by placing a flat glass under the water ||24|| . The validity 



of gaussian approximation of N(z) (BLS model) was examined by various groups and for 
a variety of systems P, ^5| , |33| . Ocko et al. have measured the reflectivity of liquid alkanes, 



over a wide range of temperatures verifying that the surface roughness is of the form given 
in Eqs. ||and |6| §3] . 



Experimentally, the reflectivity signal at each Q z from a rough interface is convoluted 
with the resolution of the spectrometer in different directions. The effect of the resolution 
along the Q z can be calculated analytically or convoluted numerically by computation. For 
simplicity, we consider that the resolution functions can be approximated as a Gaussian 
with a width of AQ Z along the Q 2 can be taken as e -< ^ 2//A< ^ with appropriate normalization 
factor ||. The resolution, AQ Z , is Q z dependent as the angles of incidence and scattering 



are varied |33|. However, if we assume that around a certain Q z the resolution is a constant 
and we measure a exp the convolution of the true reflectivity with the resolution function 
yields the following relation, 

+ AQ 2 Z (67) 



a exp a Iff 



from which the effective roughness can be extracted as follows, 



a ~ (68) 

" 'i - °LAQl 



Thus, if the resolution is infinitely good i.e., a AQ Z = the measured and effective roughness 
are the same. However, as the resolution is relaxed, the measured roughness gets smaller 
than the effective roughness. The effect of the resolution on the determination of true surface 
roughness was discussed rigorously by Braslau et al. ||. 

Diffuse scattering from liquid surfaces is practically inevitable, due to the presence of 
capillary waves. Calculation of the scattering from disordered interfaces of various char- 
acteristics were treated in |19|. In the Born approximation, true specular scattering from 



liquid surfaces exist only by virtue of the finite cutoff-length of the mean-square height fluc- 
tuations. In other words, the fluctuations due to capillary waves diverge logarithmically and 
only due to the finite instrumental resolution that true specular reflectivity is observed. The 
theory for the diffuse scattering from fractal surfaces and other rough surfaces was developed 
in PI. 
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B. Langmuir Monolayers 



A Langmuir monolayer (LM) is a monomolecular amphiphilic film spread at the air-water 
interface. Each amphiphilic molecule consist of a polar head group (hydrophilic moiety) 
and a nonpolar tail typically hydrocarbon (hydrophobic) chains |17| , |50| . Typical examples 
are fatty acids, lipids, alcohols and others. The length of the hydrocarbon chain can be 
varied chemically, affecting the hydrophobic character of the molecule. Whereas, the head- 
group can be ionic, dipolar, or with a certain shape that might attract specific compounds 
present in the aqueous solution. One important motivation to studying LM's is their close 
relationship to biological systems. Membranes of all living cells and organelles within cells 
consist of a lipid bilayers interpenetrated with specific proteins, alcohols, and other organic 
compounds that combine to give functional macromolecules that determine transport of 
matter and energy through them. It is well known that biological functions are structural, 
and structures can be determined by XR and GID. In addition, delicate surface chemistry 
can be carried out at the head-group interface with molecules from the aqueous solution. 
From the physics point of view, the LM belongs to an important class of quasi- 2D system 
with which statistical models that depend on the dimension of the system can be examined. 

Herein, results from a simple lipid, dihexadecyl hydrogen phosphate (DHDP), consisting 
of a phosphate head group (POj ) and two hydrocarbon chains attached to it, are presented. 
Figure [l(](a) displays the normalized reflectivity of a DHDP monolayer at the air-water 
interface at a lateral pressure of 40 mN/m. The corresponding electron density profile is 
shown in the inset as a solid line. The profile in the absence of surface roughness (a = 0) is 
displayed as a dashed line. The bulk water subphase corresponds to z < 0, the phosphate 
headgroup region is at > z > 3.4A and the hydrocarbon tails are at the 3.4A > z > 23.1 
A region. As a first stage analysis of the reflectivity a model SLD with minimum number of 
boxes, i = 1, 2, 3..., is constructed. Each box is characterized by a thickness di and an electron 
densities N e> i, and one surface roughness , a for all interfaces. Refinement of the reflectivity 
with Eq. [30] shows that the two box model is sufficient. In order to improve the analysis we 
can take advantage of information we know of the monolayer i.e., the constituents used and 
the molecular area determined from the lateral-pressure versus molecular area isotherm. If 
the monolayer is homogeneous and not necessarily ordered, we can assume an average area 
per molecule at the interface A, and calculate the electron density of the tail region as follows 

Ptaii = N eMil r /(Ad tai i) (69) 

where N e ta j; is the number of electrons in the hydrocarbon tail and dt a u is the length of 
the tail in the monolayer. The gain in this description is two fold; first, the number of 
independent parameters can be reduced, and constraints on the total number of electrons 
can be introduced. However, in this case, the simple relation phead = N etPhosphate / (Adh ea d) 
is not satisfactory and in order to get a reasonable fit additional electrons are necessary in 
the head-group region. These additional electrons can be associated with water molecules 
that interpenetrate the head group region which is not densely packed. The cross section of 
the phosphate head group is smaller than the area occupied by the two hydrocarbon tails 
allowing for water molecules to penetrate the head group region. We therefore introduce an 
extra parameter N# 2 o, the number of water molecules with ten electrons each. The electron 
density of the head group region is given by, 
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Phead 

(N 

e, phosphate 

+ ION H 2 o)/(Ad head ). 



(70) 



This approach gives a physical insight into the chemical constituents at the interface. In 
modeling the reflectivity with the above assumptions we can either apply volume constrains 
or equivalently examine the consistency of the model with the literature values of closely 
packed moieties. In this case the following volume constraint can be applied, 

Vheadgroup = Adhead = ^HiO^H 2 + V phosphate (71) 

where Vh 2 o ~ 30 A 2 is known from the density of water. The value of Vphosphate determined 
from the refinement, should be consistent within error with known values extracted from 
crystal structures of salt phosphate [0 . 



Another parameter that can be deduced from the analysis is, the average tilt angle, t, of 
the tails with respect to the surface from the relation, 

dtail I 'kail = COSt (72) 

where l ta u is the full length of the extended alkyl chain evaluated from the crystal data for 



alkanes [18|. Such a relation is valid under the condition that the electron density of the 



tails when tilted is about the same as that of closely packed hydrocarbon chains in a crystals 

° 3 

Ptaii ~ 0.32e/A r as observed p3| . Such a tilt of the hydrocarbon tails would lead to an 
average increase in the molecular area compared to the cross section of the hydrocarbon 
tails (A)), 

A /A = cost. (73) 

Gregory et al. found that at lateral pressure, tt = AOmN/m the average tilt angle is very 
close to zero (~ 7 ± 7°) and extract an Aq « 40. 7 A compared with a value of 39. 8A 
for closely packed crystalline hydrocarbon chains. The small discrepancy was attributed to 
defects at domain boundaries. 

The GID for the same monolayer is shown in Fig. |10|(b) where a lowest order Bragg 
reflection at 1.516A -1 is observed. This reflection corresponds to the hexagonal ordering of 
the individual hydrocarbon chains |22|,|23| with lattice constant d = 4.1144A, and molecular 



area per chain A c h ain = 19.83A 2 . Note that in DHDP the phosphate group is anchored to 
a pair of hydrocarbon chains with molecular area A = 39.66A 2 , and it is surprising that 
ordering of the head group with a larger unit cell (twice that of the hydrocarbon unit cell) is 
not observed as is evidenced from Fig. |T0"|(b). Also shown in the inset of Fig. [K](b) is a rod 
scan of the Bragg reflection. To model the rod scan in terms of tilted chains the procedure 



developed in |23[ is followed. The structure factor of the chain can be expressed as 



F cha UQ'^Q' z ) = F ^ f^jw ( 74 ) 

where F{Q' 1 ) is the in-plane Fourier transform of the cross section of the electron density of 
chain weighted with the atomic form factors of the constituents. The second term accounts 
for the length of the chain and is basically a Fourier transform of a one dimensional aperture 
of length I. If the chains are tilted with respect to the surface normal (in the y-z plane) by 
an angle t, the Q' should be rotated as follows, 
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Q' x = Q x cost + Q z sint 

Q y — Qy 

Q' z = —Q x sin t + Q z cost (75) 

For small chain-tilt angles this rotation mainly affects the Q z dependent part of the chain 
structure factor, since the Qx_ changes. Applying this transformation to molecular structure 



factor, Eq. |74|, and averaging over all six domains (see more details in Ref. with the 
appropriate weights to each tilt direction we find that at 40 mN/m the hydrocarbon chains 
are practically normal to the surface consistent with the analysis of the reflectivity. 

In recent Brewster Angle Microscopy (BAM) and X-ray studies of Cgo-propylamine 
spread at the air- water interface (see more details on fullerene films ||55||), a broad in-plane 
GID signal was observed |l6j. The GID signal was analyzed in terms of a 2D radial dis- 
tribution function that implied short range positional correlations extending to only few 
molecular distances. It was demonstrated that the local packing of molecules on water is 
hexagonal, forming a 2D amorphous solid. This is a detailed study demonstrating how to 
analyze homogeneously disordered 2D system by combining X-ray scattering techniques and 
visible light microscopy. 



C. Surface Crystallization of Liquid Alkanes 



Normal alkanes are linear hydrocarbon chains (CH 2 ) n terminating with CH 3 groups 
similar to fatty acids and lipids that by contrast posses a hydrophilic head groups at one 
end. Recent extensive X-ray studies of pure and mixed liquid alkanes |58|-|60]|34| reveal a rich 
and remarkable properties near their melting temperature, Tf. In particular, a single crystal 
monolayer is formed at the surface of an isotropic liquid bulk up to ~ 3°C above Tf for a 
range of hydrocarbon number n. The surface freezing phenomena exists for a wide range 
of chain lengths 16 > n > 50. The molecules in the ordered layer are hexagonally packed 
and show three distinct ordered phases: two rotator phases, one with the molecules oriented 
vertically (16 > n > 30) and the other tilted toward nearest neighbors. (30 > n > 44). The 
third phase (44 > n) orders with the molecules tilted towards next-nearest neighbors. In 
addition to the 2D Bragg reflections observed in the GID studies, reflectivity curves from 
the same monolayers were found to be consistent with a one box model of densely packed 
hydrocarbon chains, and a thickness that corresponds to slightly tilted chains. This is an 
excellent demonstration where no other technique but the X-ray experiments carried out at 
a synchrotron could be applied to get the detailed structure of the monolayers. Neutron 
scattering from this system would have yielded similar information, however the intensities 
available today from reactors and spallation sources are smaller by at least a factor of 10 5 
counts/sec for similar resolutions and will not allow observation of any GID signals above 
background levels. 



D. Liquid Metals 



Liquid metals unlike dielectric liquids consist of the classical ionic liquid and quantum 
free electron gas. Scattering of conduction electrons at a step like potential (representing the 



22 



metal- vacuum interface), give rise to quantum interference effects and lead to oscillations of 



the electron density across the interface |26j . This effect is similar to the Friedel oscillations 
in the screened potential arising from the scattering of conduction electrons by an isolated 
charge in a metal. By virtue of their mobility, the ions in a liquid metal can in turn 
rearrange and conform to these oscillations to form layers at the interface, not necessarily 
commensurate with the conduction electron density [[|l|]. Such theoretical predictions of 
atomic layering at surfaces of liquid metals were known for a long time and were only 
recently confirmed by X-ray reflectivity studies for liquid gallium and liquid mercury |3iy40| . 
X-ray reflectivities of these liquids were extended to Q z ~ 3A showing a single peak that 
indicates layering with spacing on the order of atomic diameters. The exponential decay for 
layer penetration into the bulk of Ga (6.5A) was found to be larger than that of Hg (~ 3 



A). Figure [11] shows a peak in the reflectivity of liquid Ga under in situ UHV oxygen free 



surface cleaning | 40| . The normalized reflectivity was fitted to a model scattering length 



density shown in Figure |TT|(b) of the following oscillating and exponentially decaying form 

p(z)/ Ps = erf[(z - zo)/a] + 9(z)A sin(27rz/rf)e" 2/§ (76) 

where 9(z) is a step function, d is the inter-layer spacing, £ the exponential decay length, 
and A an amplitude. Fits to this model are shown in Fig. |y] with d = 2.56A, £ = 5. 8 A. The 
layering phenomena in Ga showed a strong temperature dependence. Although liquid Hg 
exhibits layering with a different decay length the reflectivity at small momentum transfers, 
Q z are significantly different than that of liquid Ga indicating fundamental differences in the 
surface structures of the two metals. The layering phenomena suggests in-plane correlations 
that might be different than those of the bulk, but had not been observed yet with GID 
studies. 
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APPENDIX A. 



p-polarized X-ray beam 



A p-polarized X-ray beam has a magnetic field component that is parallel to the stratified 
medium (along the x-axis, see Fig. p, and straightforward derivation of the wave equation 
|7] yields 




kl - V(z) 



B = 0. 



(77) 
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By introducing a dilation variable Z such that, dZ = edz Eq. [77| for B can be transformed 
to a form similar to Eq. 12 



+ | fizZM 1B .O. ,78, 

The solution of Eq. for an ideally flat interface in terms of r p (k z , s ) and t p (k Zta ) is then 
simply given by 

^ = tiiti' tp(h " s> = w <79) 

The critical momentum transfer for total external reflectivity of the p-type X-ray beam is 
Q c = 2k c = 4^71 p s , identical to the one derived for the s-type wave. Also, for 2kf 3> Q z 3> 
Q c , (kf is defined below), Rf(Q z ) can be approximated as 

MQ.)~(&)U^-) a - (so) 



2Q Z J VI + e 

The factor on the right hand side is one for all practical liquids, and thus the Born approxi- 
mation is basically the same as for the s-polarized X-ray beam (Eq. ^7|) . The main difference 
between the s-type and p-type waves occurs at larger angles near a Brewster angle that is 
given by 9b = sin _1 (/cf /ko). At this angle, total transmission of the p-type wave occurs 



r p (k ZjS ) = 0). Using Eqs. Ill and 179], kj can be derived 



k~ 1 , * 

-r= i si 

kn JO. — A-rrn / h? 



o yj2 - 4np a /k% 

The Brewster angle for X-rays is then given by 6b = sin _1 (/cf /ko) ~ 7r/4. This derivation 
is valid for solid surfaces, including crystals, where the total transmission effect of the p- 
polarized wave at a Bragg reflection is used to produce polarized and monochromatic X-ray 
beam. 
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FIGURES 



FIG. 1. The geometry of incident and scattered beam in specular reflectivity (a) and in a 
non-specular scattering(b) experiments. 

FIG. 2. Calculated reflectivity curves for external (solid line) and internal (dashed line) scatter- 
ing from an ideally flat interface versus momentum transfer given in units of the critical momentum 
transfer, Q c = 4(-7rp s ) 1//2 . The dotted line is kinematical approximation (Qc/^Qz) 4 - The lower panel 
shows the amplitude of the wave in the medium for external (solid line) and external reflection 
(dashed line). 

FIG. 3. Calculated reflectivities from H2O and liquid mercury (Hg) showing the effects of 
absorption and surface roughness. The absorption modifies the reflectivity near the critical mo- 
mentum transfer for mercury with insignificant effect on the reflectivity from H2O. The dashed line 
shows the calculated reflectivity from the same interfaces with root mean square surface roughness, 
a = 3A. 

FIG. 4. An illustration of a continuous scattering length density, sliced into a histogram. 

FIG. 5. Calculated reflectivities for two films with identical thicknesses but with two distinct 
normalized electron densities pi (solid line) and 1 — pi (dashed line) and corresponding calculated 
reflectivities using the dynamical approach (a = 0). The two reflectivities are almost identical 
except for a minute difference near the first minimum (see arrow in figure) . The Born approximation 
(dotted line) for the two models yields identical reflectivities. The inset shows the normalized 
reflectivities near the first minimum. As Q z is increased the three curves converge. This is the 
simplest demonstration of the phase problem, the non-uniqueness of models where two different 
potentials give the same reflectivities. 

FIG. 6. (a) Illustration of wave paths for exterior (a) and interior(b) scatterer near a step-like 
interface. In both cases the scattering is enhanced by the transmission function when the angle of 
the incidence is varied around the critical angle. However, due to the asymmetry between external 
and internal reflectivity the rod scan of the final beam modulates the scattering differently as is 
shown on the right hand side in each case. 

FIG. 7. Monochromator geometry to tilt a Bragg reflected beam from the horizon on a liquid 
surface. Two possible tilting configurations about the primary beam axis and about an axis along 
the surface of the reflecting planes are shown. 
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FIG. 8. Superposition of the reflected-beam (circles) below the critical angle and direct beam 
(triangles), demonstrating total reflectivity of X-rays from the surface of water. Severe surface 
roughness reduces the intensity and widens the reflected signal. A reduction from total reflectivity 
can also occur if the slits of the incident beam are too wide, so that the beam-footprint is larger 
than the surface sample. 

FIG. 9. Experimental reflectivity from the surface of water. The dashed line is the calculated 
Fresnel reflectivity from an ideally flat water-interface, Rf- The Normalized reflectivity versus Q 2 Z 
is fitted to a the form R/Rp = e~(Q^ a ) , demonstrating the validity of the capillary-wave model 
[Buff, et al., 1965]. 

FIG. 10. (a) Normalized X-ray reflectivity from Dihexadecyl-phosphate (DHDP) monolayer 
at the air-water interface with best fit electron density, N e , shown with solid line in the inset. 
The calculated reflectivity from the best model is shown with a solid line. The dashed line in the 
inset shows the box model with no roughness a = 0. (b) A diffraction from the same monolayer 
showing a prominent 2D Bragg reflection corresponding to the hexagonal ordering of individual 
hydrocarbon chains at Q B _L= 1.516 A" 1 . The inset shows a rod scan from the quasi-2D Bragg 
reflection at Q±, with a calculated model for tilted chains denoted by solid line (see text for more 
details). 

FIG. 11. Upper panel shows measured reflectivity for liquid Ga. Data marked with X were 
collected prior to sample cleaning whereas the other symbols correspond to clean surfaces (for 
details see Regan et al., 1995). Calculated Fresnel reflectivity from liquid Ga surface convoluted 
with a surface roughness due to capillary waves (a = 0.82A), and the atomic form factor for Ga 
is denoted with a solid line. The lower panel shows the normalized reflectivity, with a solid line 
that was calculated with the best fit by an exponentially decaying sine model shown in the inset 
(Courtesy of Regan et al. 1995 
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TABLES 



TABLE I. Electron number density, SLD, critical angles and momentum transfers, and absorp- 
tion term for water and liquid mercury. 





N e (e/A 3 ) 


Ps (A 2 x io- 5 ) 


Q C {A- 1 ) 


a c (deg.) 


/3(xl0- 8 ) 










for A = 1.5404 A 




H 2 


0.334 


0.942 


0.02176 


0.153 


1.2 


Hg 


3.265 


9.208 


0.06803 


0.478 


360.9 
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